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Abstract. - This letter provides the first analytical estimation of the probability density function 
(PDF) tail of the interface width in the Kardar-Parisi-Zhang equation in the regime of strongly 
non-linear growth. We find that the PDF tail is proportional to exp{— cuij''^}. In addition, the 
effect of spatial dimensions on the PDF tail scaling is dicussed. The PDF tail is computed by 
the instanton method within the Martin-Rose-Siggia framework using a careful treatment of the 
non-linear term. This gives a novel approach to understand the rightmost PDF tail of the interface 
width distribution and the analysis suggests that there is no limit in the upper critical dimension. 



Introduction. — In nature there are niany important 
phenomena that are driven far from equilibrium by insta- 
bilities or by external forces. Examples are diverse from 
forest fires to interstellar turbulence, which is constantly 
stirred by supernova explosions. A proper description and 
understanding of the multiscale interactions that are re- 
sponsible for the inevitably complex dynamics in these 
nonequilibrium systems remains a significant challenge in 
classical physics. 

The out of equilibrium interfacial growth is another ex- 
ample that has attracted much attention during recent 
years. A description of these growth processes that have 
been widely recognized is a Langevin like equation by 
Kardar-Parisi-Zhang (KPZ) [T]. The KPZ equation is one 
of the simplest non-linear generalizations of the diffusion 
equation and is thus connected to many other areas of non- 
equilibrium dynamics such as the Burgers turbulence [2]- 
[3] , driven diffusion and dissipative transport |4l as well as 
flame front propagation 0. 

The KPZ equation has been studied extensively, how- 
ever there are some remaining controversial issues, in par- 
ticular estimates of the upper critical dimension are in the 
range dc = 2.8 — oo [B]- [13]. The reason for searching for a 
particular value of the upper critical dimension (i.e. below 
dc there is no phase transition) is the hope of systematic 
expansions in dc — d in analogy with equilibrium critical 
phenomena. 

The purpose of the present work is to provide a sta- 



tistical theory of interfacial growth in the strongly non- 
linear growth phase and thereby shed light on the elu- 
sive finite upper critical dimension in the KPZ equation. 
We compute the tail of the probability density function 
(PDF) of the interface width using the instanton method 
in the Martin-Siggia-Rose framework [21]. The instanton 
method is a non-perturbative way of computing the PDF 
tails [22]- [17]. The PDF tails can be viewed as a transi- 
tion amplitude from a quiescent state (where no growth 
occurs) to a final state determined by a coherent struc- 
ture, from which the PDF tails are computed through a 
path-integral. Similarly to the Burgers equation the KPZ 
equation supports a coherent structure solution that can 
be used in the calculation as the path with highest prob- 
ability and subsequently the path-integral can be solved 
using the saddle-point method. The instanton predictions 
of the PDF tails have been shown to agree very well with 
numerical simulations in a model of self-organization of 
sheared fiows [25] and heat fiux |26j. 

The PDF has been computed in a similar manner pre- 
viously in Refs. [TJ]- [TB], however all these results relies 
on the assumption of a weakly non-linear system. In the 
present setting we focus on the effects in the intermittent 
or strongly non-linear regime. 

We find that the PDF exhibit heavier tails than a Gaus- 
sian distribution whereas the tails are subdominant to 
the exponential distribution. Under the assumption of 
isotropic growth in d-dimensions we find a smooth vari- 
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ation as the dimensionality d is increased. This suggests 
that there exists no upper critical dimension and this is in 
accordance with the numerical work in Ref . [6] . 

The letter is organized as follows. First, the Kardar- 
Parisi-Zhang model is introduced and the PDF tail is com- 
puted using the instanton method. Then the result of the 
instanton PDF tails are corroborated in the one dimen- 
sional case through a derivation using the Fokker-Planck 
method. Finally the generalization to higher spatial di- 
mensions are introduced and the letter is concluded by 
results and discussion. 

The PDF tail of the KPZ equation. Lateral 
interface growth can be described by the Kardar-Parisi- 
Zhang (KPZ) model [I], that is the time evolution of the 
height h is proportional to the square of the height gradi- 
ent modified by diffusion and a stochastic forcing, 



at 



(1) 



Here / is a white noise forcing with a short correlation 
time modeled by the delta function as 



1 



:6(t-t')cxp{-{{x-x'/ey} 
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5{t - t')K,^{x - x'), 



(2) 



and (/) = 0. The angular brackets denote the average 
over the statistics of the forcing /. The 1-1-1 dimensional 
KPZ equation is equivalent to the noisy Burgers equation 
by the relations 



udx. 



(3) 
(4) 



Using the relations in Eq. ([3])-(|4]), we now find the noisy 
Burgers equation to be 



— + 2AuVu - j/V^u = -V /. 
dt 



(5) 



The Burgers equation is known to support the ramp and 
shock- like stationary coherent solutions [5] . We will adopt 
the ramp solution of the form u on x, using the relations 
between Burgers and KPZ we can find the solution for h 
of the form h — ax^ + b where the time evolution of h and 
the non- linear term both will be proportional to x^. The 
constants a and b can be determined by inputting the trial 
function into KPZ yielding a = and b = j^. Using the 
parabola solution as a coherent structure for the instanton 
method we will now compute the PDF tails. 

We compute the PDF tails of the mean square height 
fluctuations (the interface width) where we give special 
attention to the effect of dimensional scaling by using the 
instanton method. Here the mean square height fluctua- 
tions are defined as 



W2 



(6) 



where is the area of the substrate with characteristic 
linear dimension L, and h = J2r ^(^lO/^i is the average 
height of the surface. The PDF tails of W2 are expressed 
in terms of a path-integral using the Gaussian statistics 
of the forcing. The optimum path is then associated with 
the creation of a short lived coherent structure (among 
all possible paths) and the action is evaluated using the 
saddle-point method on the effective action. The saddle- 
point solution of the dynamical variable h{x,t) of the form 
h{x, t) = F{t)(j){x) is called an instanton if F{t) — Q aXt = 
—oo and F{t) 7^ at t = 0. Note that, the function (j){x) 
here represents the spatial form of the coherent structure. 
We will first consider the d — 1 case and then generalize 
the found PDF to arbitrary dimension d. The probability 
density function of W2 can be defined as 

P{W2) = {SiMih) - W2)) = J d^e'^^U^, (7) 



where 



k = (exp(-zCM(/i))). (8) 



Here M{h) is the general expression for the m-th moment 
(/i™), however we will restrict this study to the second 
moment. Following Ref. [52] the integrand can then be 
rewritten in the form of a path-integral as 



(9) 



where effective action S'^ of the KPZ equation is expressed 
as 



I dxdth [ ^ - \[Vhf - vV'^h 



+ J dxdx'dth{x,t)K^{x ~ x')h{x' ,t) 

+ j dxdth^dit). (10) 

Utilizing the instanton function h{x,t) = F{t)4>{x) the 
action can now be recast into, 

= -i I d^i{ciF - c^XF^ - C2vF^ 



(11) 



C4 j dtfi^ 

J dtcsF^Sit). 



The conjugate variable is denoted h — fi{t)(j) and we have 
used the following definitions of the constants. 



Co 



ci = 



C2 



dx(l){x)(W<l>)^{x) 
dx(l){x)4>{x), 
dx(f>{x)'V'^ (p{x) , 



(12) 
(13) 
(14) 
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C3 



C4 



dx(j)^{x), 
dxdy(f>{x)(t>{y), 



(15) 
(16) 



Note that the constant C4 is evaluated for small values 
of e in Eq. ^ and in the higher dimensional case C2 will 
change with the dimension d and C2 ^ co . We compute the 
first variational derivatives to minimize with respect to 
F and /i in order to find the path with highest probability 
identified by the instanton or the extremum of the action 
as, 

= (ci^ - coAF^ - jycaF) + 2c4fi = 0, (17) 

s s 

= (-ci/i - 2coXFn - i'C2fi) + 2i£,c3FS{t) = 0. 

O-C 

(18) 

Eqs. (fT7|) - ([T8| constitutes a dynamical system for the 
instanton time function F and its conjugate /i. We proceed 
by solving these equations for t < and matching the 
solution at t = 0. Note that the instanton solution F 
rapidly grows at t = with increasing ^ while it vanishes 
as i — —00. We start by computing an additional relation 
for the time evolution of the conjugate variable fj, expressed 
in the real variable F by taking the time derivative on Eq. 

(ED, 



CiF - 2coXFF - VC2F 



-2zc4/i. 



(19) 



We now substitute and fi in Eq. (|19|) by using Eqs. 
lfT7|) and (|18p yielding a second order non-linear differen- 
tial equation for _F, 



ciF = v^c^^ + 3i/coC2AF^ + 2cl\^F^ 



(20) 



By setting v — F we find that we can rewrite the time 
derivative as v{dv/dF) and we can now perform the first 
integration, 

cfv'^ = ly^cjF^ + 2ucoC2\F^ + c^A^F^ 

= F^{vc2 + cq\FY. (21) 

Below we will utilize the relation between F and F that 
can be expressed as. 



ciF = ±F(j/C2-f coAF). 



(22) 



We can now easily determine the instanton time depen- 
dence from the separable differential Eq. (|22|) . 



m = 



VC2 



H - 

VC2 



CqA' 

+ coAF(O) 



F(0) 



(23) 
(24) 



Note that the equation for F[t) gives F(0) at t = 0. We 
now have to determine a value of F at t = as a function 



of ^. Thus we integrate Eq. (IT^ over (— e, e) and use Eqs. 
([T7)) and (HH) with the negative sign to obtain. 



2^^f. 

CqCiA 



(25) 



The path-integral will now be computed using the saddle- 
point method in the limit of ^ — >■ cxi. First we have to 
evaluate the ^-dependence of the action S^. In the limit 
of ^ — > 00, becomes 



Si 



dt[i [cxF - coXF^ - vc2F^ 
C4 / dt^^ 



iS. j dtc3F^5{t) 

dt{ciF - coXF^ 



1 

4C4 
C4 Jo 



vc2F)^ +i^C3F^{0) 



F(0) 



Cl 

— \lyc2 

C4 



dF{vc2F + coXF^) + i^C3F^{0) 
F2(0) ^ ^ ,F3(0)\ 



3 CqC^A^ 

4. 44 



CqA- 
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3 y 



iCc3F'{0) 



3 CqC^A^ 



Now let 



4 c^c^ 



3 CqC^A^ 



(26) 



(27) 



and the action becomes = — iCC^- The tail of the PDF 
is then found by performing the ^-integral in Eq. ([7]) by 
the saddle point method in the limit W2 00. It is later 
shown that this corresponds to ^ — > cx). 



Piw2) 



d^e 



(28) 
(29) 



We evaluate the ^-integral using the extreme point 
/'(Co) = and ^2 = -i?/(3C) of /(^ = *C^2 + ^C?^ for 
the saddle-point method. This results in the PDF of W2 
as, 



Piw2) 



O 3/2 

exp-^ =■ 

^ 3V3C 



} 



(30) 



where ^ is determined by Eq. ((27l 



Computing the PDF tails using the Fokker- 
Planck equation. — In order to validate our found PDF 
we perform a similar analysis for d = I where we compute 
the PDF using the Fokker-Planck (FP) method [22], [27]. 
However, it seems to be very difficult to generalize the 
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PDFs found using the FP method to arbitrary spatial di- 
mensions. Furthermore, the PDF is found to have the 
same fundamental exponential form whereas the coeffi- 
cients differ. We assume that we can write the functions 
{h and /) in Eq. ([T|) as 



h{x,t) = (p{x)F{t), 
f{x,t) = 



(31) 
(32) 



with (p as before and {g{t)g{ti)) — GS{t—ti). By substitut- 
ing this into Eq. (1) we find (neglecting the dissipation), 



F = AXaF^ + g{t). 



(33) 



Furthermore we write the generating function (z) and the 
PDF (P) as 



(34) 



(z) = j dFP{F)e^^^ ^P{il (35) 
P(F) ^ I die'^^{z). (36) 



Here the (•) is the mean value integral over the forcing /. 
To compute the PDF from the FP equation we have to 
determine the time evolution of the generating function z 
using Eq. ([M)) . 



dz 



i^Fe'^^ 

i^iiXaF'^ + g{t))z. 



(37) 
(38) 



The time evolution of P can now be found using Eq. 
- dSni) as 

dt dF^ ' dF^ 



(39) 



In the case of stationary PDF the differential equation is 
separable and we find 

PiF) ^ Poe-"^^", (40) 

where this can be rewritten in terms of and we find 

(41) 



^/^On ^ 4Aa /r,2>3/2 



This is the same stretched exponential dependence as was 
found using the instanton method whereas the coefficient 
is different in the case of one spatial dimension. 

Dimensional scaling of the PDF tail. In this 
letter we will model the effects of several spatial dimen- 
sions for the PDFs found from the instanton method as 
being isotropic in all directions. Note that the KPZ has 
two terms involving the gradient operator where we as- 
sume V4>{xi,X2, Xd) ~ k4>{xi,X2, Xd) where kd is the 
spatial length scale and (f) is the spatial coherent structure 
in d dimensions. In the isotropic case the terms become 
V^</) = dk'^4> and (V(^)^ = dk'^cj)'^. The generaUzation to 



the anisotropic case is straightforward, however, it does 
not add significantly to the understanding of scalings of 
dimensionality in the PDF. The coefficients in the found 
PDF transforms as. 



cod = dkjco, 



C2d 



Cld = Ci, 

dkdC2- 



(42) 
(43) 
(44) 



Here the coefficients cod,cid,C2d represents the trans- 
formed coefficients in the d-dimensional case. We will now 
show some results of the predicted PDF tails combining 
Eq. ^ with Eqs. (|42l) -(|i4l). 

We will now elucidate on the effects of in particular 
dimensionality but also on the significance of the non- 
Gaussian tails by computing the dimensional dependence 
of the PDF tails. For the case of dimension = 1 we will 
compare the analytical prediction in Eq. (j30p with an ex- 
ponential tail, P{w2) ^ e^'^'^'^^, which we motivate below. 
The expression for the PDF tails can be written as 



P{W2) 



cdw\ 



3/2 



(45) 



where c is a numerical coefficient and d is the spatial 
dimension. Note that the power 3/2 in the exponential 
comes directly from the non-linear interaction term in the 
KPZ equation whereas the factor cd arises from the linear 
part of the equation. This suggests that KPZ supports 
harmonic growth in all spatial directions and the physical 
dimensions are weakly correlated with a smooth behaviour 
of the scaling function as the dimensionality changes. Fur- 
thermore, it is important to note that to have a finite up- 
per critical dimension, the scale lengths kd in Eqs. 
and have to change. 

Much of the knowledge about the KPZ equation come 
from simulations on systems, typically lattice models, that 
belong to the KPZ universality class. The problem with 
such lattice models is that given a specific model, such as 
the restricted solid-on-solid model [5S], the KPZ parame- 
ters, A and v are fixed and the non-linear parameter. A, is 
often quite small [29] , which result in PDF tails that are in- 
distinguishable from exponentially decaying PDF tails [IB] 
or result in stretched exponentials with exponent slightly 
smaller than one [301 ■ This is is contrast to our exponent of 
3/2, which holds in the strongly non-linear regime. Due to 
numerical instabilities, which seem to grow with increasing 
A, it is a very challenging task to solve the KPZ equation 
by numerical integration. Thus, we are not aware of any 
numerical estimations of the solutions to the KPZ equa- 
tion in the strongly non-linear regime with which we could 
compare our analytical results. 

Results and discussions. — In finding the analytical 
solution for the PDF we assumed that we solved for values 
in the rightmost part of the tail due to taking ^ — >■ c» and 
using the saddle-point method whereas here we will show 
the PDF for all values of W2/{w2) > 1. The constant 
c in Eq. (j45|) includes an unknown forcing strength and 
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Fig. 1: (Color online). The scaling function of the width dis- 
tribution P{w2) as a function of W2/{w2) (linear scale) where 
the results are shown for exponential PDF (red) and instanton 
PDF (black). 



Fig. 2: (Color online). The scaling function of the width dis- 
tribution W2P{w2) as a function of W2/iw2) (log-linear scale) 
where the results are displayed for exponential PDF (red) and 
instanton PDF (black). 



can be estimated from simulations, e.g. those reported in 
Ref. [6] . In Figure 1 we show the scaling function of the 
width distribution P{w2) as a function of ^2/(^2) where 
quite small differences are present. Note that the PDF 
is normalized in such a way that the total probability is 
unity and that -P(l) ~ 1 whereas the values of 102/(^2) are 
arbitrary. The maximum value of the PDF is determined 
by the strength of the forcing. Figure 2 displays again the 
scaling function P{'W2) as a function of 102/ {W2) however 
this time in lin-log scale which highlights the significant 
differences for small probabilities. 

To this end, in this paper we have computed the PDF 
tails in the KPZ model and evaluated the dependence 
of spatial dimensionality. The PDF tails of the form 
~ exp{— 0^2^^} (where c is a numerical coefficient de- 
pendent on the dimension d) are relevant for intermittent 
events in the strongly non-linear growth regime. More- 
over, the PDF tails are significantly alleviated compared 
to a Gaussian distribution whereas the tails are subdomi- 
nant to the exponential distribution. Of particular inter- 
est are the effects of higher dimensions on the coefficient 
c. Here it is shown that using an isotropic growth model 
the coefficient c changes smoothly as d increases and thus 
suggests that there is no upper critical dimension dc- 

^ ^ ^ 
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